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Abstract 

We map self-avoiding random walks with a chemical potential for writhe 
to the three-dimensional complex 0{N) Chern-Simons theory as — > 0. We 
argue that at the Wilson-Fisher fixed point which characterizes normal self- 
avoiding walks (with radius of gyration exponent u ~ 0.588) a small chemical 
potential for writhe is irrelevant and the Chern-Simons field does not modify 
the monomer- monomer correlation function. For a large chemical potential 
the polymer collapses. 
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1 Introduction and Summary 



Aside from their mathematical interest, self-avoiding random walks are good models for 
polymers which are much longer than their persistence length. In this paper, we study 
how these self-avoiding random walk models can be adapted to study the effects of a 
topological constraint on the polymer. Recent progress in experimental technique has, 
for the first time, allowed the direct, real-space observation of polymer conformations ]![. 
These studies, as well as experiments in which polymers have been subjected to torsional 
constraints [Q, have focussed primarily on long biopolymers. While self-avoidance is 
important in long biomolecules such as actin or DNA, these molecules, in addition, have 
a twist rigidity. If the twist were decoupled from the conformational changes of the 
polymer then twist rigidity would be unimportant. However, for a closed polymer, such 
as a plasmid DNA loop, there is a topological invariant, namely the linking number of 
the two DNA backbones. Fuller's now ubiquitous [^, ||, 0] relation Lk = Tw + Wr 
provides a connection between the topological linking number Lk, the total amount of 
twist in the polymer Tw and the writhe of the curve Wr — a scalar which captures some 
geometrical information about the conformational deformations of the loop. The way 
that topological constraints can affect the conformational statistics of polymers is now 
accessible to experimental investigation. Studies of DNA force-extension curves with 
constrained Lk have discovered a coupling between an excess in linking number and the 
backbone configuration 0]. The question arises: could there be a new universality class 
describing the scaling regime of topologically constrained polymers? 

There is currently a wealth of information on the statistics of self-avoiding random 
walks. In particular, a most remarkable mapping by de Gennes between the critical 
behavior of the scalar 0{N) model as — and the behavior of self-avoiding walks 
has allowed, via any number of schemes, a method of calculating the radius of gyration 
exponent u which measures how the radius of gyration scales with the polymer length: 
Rg ^ L^ . It is interesting then to consider the effect of the topological constraint on the 
radius of gyration of the closed polymer loop. The topological constraint can be added 
to the standard analysis via a gauge field, as noted by Brereton and Shah 0. 

As we will argue, polymers with a small chemical potential for Lk fall into the same 
universality class as the usual self-avoiding random walk. The experimentally observed 
coupling arises from the non-universal details of the polymer and not from any general 
long-wavelength properties. A chemical potential for Lk only becomes noticeable when 
it exceeds a certain threshold, at which point the chain collapses 0]. 

So far we have focussed our attention on a chemical potential for Lk. In what follows, 
we turn our attention to the writhe of a curve. The introduction of a Chern-Simons 
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gauge field allows us to naturally and easily incorporate a chemical potential for writhe 
1^. In section 2 we will show that a chemical potential for writhe is equivalent to a 



chemical potential for link. What is more, our model will only have one backbone to 
keep track of. 

In section 3 we map the statistics of self-avoiding walks with a chemical potential 
for their writhe onto the N complex scalar Chern-Simons theory. We then argue 



in section 4 that at the self-avoiding fixed point controlled by z/ ~ 0.588 |Tl| the effect 
of writhe is irrelevant for a small chemical potential. From this we conclude that there 
is no coupling between the conformations of a polymer and its writhe in the scaling 
regime. Finally, in section 5 we calculate moments of the writhe and find that the field- 
theoretic analysis disagrees with some rigorous results ll2|, We understand this to be 
a consequence of the writhe "hiding" from the long wavelength field theory in persistence 
length size, plectonemic regions separated from each other. 



2 Integrating Out the Twist 

The theory we develop in section 3 incorporates a chemical potential for the writhe of a 
closed polymer loop. While it would be difficult, experimentally, to add such a chemical 
potential, a chemical potential for linking number could easily be arranged. For example, 
DNA plasmids in the presence of topoisomerase and ATP will tend to have some non- 
zero excess linking number due to the bias of the added enzyme. We show in this section 
that a chemical potential for writhe is equivalent to a chemical potential for the linking 
number. 

For simplicity we assume that the long-wavelength behavior of the polymer is 
adequately described by a simple continuum model. The free energy for this model 
is determined by a bending stiffness n and a twist rigidity C[^: 




where R(s) describes the conformation of the chain backbone, fl{s) is the twist degree 
of freedom and s is the curves arclength. The partition function for this model with a 
chemical potential for the linking number is then given by 

Z = J [dR] [dQ] exp |-F[R, Q] - g^ (^Wr[R] + ^ds^ | , (2.2) 

where the energy is measured in units of kj^T. We have used Fuller's formula Lk = 
Tw + Wr explicitly to divide the linking number into the integral of the twist density Q 
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and the writhe, where the writhe of the curve is 



^ ^ ^nJ J \ ds ds' J |R(s) -R(sO|' 

The twist degree of freedom can be integrated out to yield a new effective partition 
function with a chemical potential for writhe: 



= I m exp / ds - g'Wr[R] 



(2.4) 



Unlike the twist, the writhe is non-local and can not be integrated out so easily. 

In other models it may be somewhat more cumbersome to integrate out the twist 
degree of freedom. As an example we can consider a ribbon-like model. Once again we 
will parameterize the backbone by R(s), but we will use a unit vector u(s) to keep track of 
the twist. We take this vector to be perpendicular to the backbone tangent, for instance, 
in DNA u could point from the central backbone to one of the two sugar-phosphate 
strands. In terms of this parameterization, the twist density is 



ds \ ds . 



(2.5) 



The partition function for this modified model is 

Z = J [dR] [dn] [du] exp |-F [R, Q] - g"^ (^Wr[R] + ^ J Qds 

6[u' -1]6 



xS 



^ dR du^^^ 

ds \ ds , 



dR 



ds 



(2.6) 



Upon integrating out u. 



[dR][dfi]det"^ 



dR 

ds 



5u 



xexp|-F[R, n]~g^ [wr\R] + ^ J dsQ^^ 



(2.7) 



The determinant is an isotropic functional of the backbone tangent. This functional can 
be expanded in powers of dR/ds. Since at distances large compared to the persistence 
length any local isotropic model flows to the random walk fixed point we will, again, only 
consider the effective partition function given in ( p.4|) . 
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3 Mapping to ^ Scalar Chern-Simons Theory 

In this section, we map the — > 0, 0{N) complex scalar field theory with a Chern- 
Simons gauge field to a self-avoiding random walk with a chemical potential for writhe. 
The Chern-Simons term is introduced because after functionally-integrating it out we 
obtain an expression for the writhe of the self-avoiding random walk described by the 
scalar theory |T^. A similar procedure to this one was used to study the area of walks 
in two-dimensions by Cardy and the linking number of a pair of walks where one 
walk fills space , though in both those cases the gauge field had the standard Maxwell 
potential F^^. 

We start by considering the Wilson loop along a closed path F for a Chern-Simons 
gauge theory with coupling g. Namely 

C[T] = ^exp|zj(rf/?^A^|^ (3.1) 
= J[dA]expi^i yrfs^-A^(R(s))|exp|-^ J d^x A-V x exp {-Sg{} , 

where we have added the gauge fixing term, 

5gf = ^/A(V-A)2. (3.2) 

Though the Wilson loop is gauge-invariant, in practice we must select a gauge to compute 
it. We choose the Landau gauge (A — ^ 0) so that the Chern-Simons propagator takes on 
a simple form: 

{A{q)A,i-q)) = + ^ 2/^. (3.3) 

As we shall show, the antisymmetry of the propagator in this gauge simplifies the 
renormalization group analysis by eliminating a number of Feynman graphs. 
Writing 

(R(s)) = J d\ 5^ [x - R(s)] A^(x) (3.4) 
we can perform the functional integration in (|3.1| ) to find 



exp 



|i£cii?^A^|^ =exp{-/Wr[F]}. (3.5) 



The writhe of any closed path is therefore determined by the line integral of the Chern- 
Simons gauge field. 

To include self-avoidance, we employ a generalization of de Gennes' mapping of self- 
avoiding random walks onto the N —' 0, 0{N) real scalar field theory It is useful 
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instead to consider the ^ 0, 0{N) complex scalar field theory which is equivalent ||15 . 



Recall that the mapping proceeds by considering the high temperature series for the N 
component complex Ising model with Hamihonian 

H =J24-Sj + c.c (3.6) 

where denotes nearest neighbors and Sj is an A^-component complex vector field 

with sJ-Sj = A^. The discrete version of the Wilson line for a path from to x is 

C[r] = (exp{tY.AA'j (3.7) 

where Aij is a field defined on the link connecting sites i and j and the sum is along 
the path F which connects and x. In order to weight each path by its writhe, we 
can include in the series expansion the corresponding Wilson line and integrate out A 
to give each path a chemical potential g"^ for writhe. Incorporation of this phase is 
straightforward — we replace the Hamiltonian in (|3.6|) with a modified one: 



H = J2 4- ^jUij + C.C., (3.8) 
where Uij = exp{iAij}. We recognize this Hamiltonian as that for a lattice gauge theory 



Iq] , where the local gauge symmetry is 



Si ^ Sic''^ (3.9) 
where 6i is arbitrary. In the continuum limit the transformation of Uij becomes 



A ^ A + 89, the usual gauge transformation. Taking the continuum limit [0 of ( p.9| ) 
and including the Chern-Simons term leads us to the 0{N) complex scalar Chern-Simons 
action 



S= £x 



(9-zA)^|'-^|V^|2 + M|^|^ + ^A-VxAj , (3.10) 

where ip is an iV-component, complex scalar field, u is the monomer-monomer repulsion 
and n is the chemical potential to add a monomer or step. Note that, though the scalar 
field theory has a U{N) symmetry, we are only gauging the diagonal U{1). There is no 
quantization condition on g as is typical in non-Abelian Chern-Simons theories fl^ or in 
the presence of magnetic monopoles [|1^, , which are not part of this theory. 
Recalling the mapping p, [1^, the correlation function is 

Gr(x)= hml(v^t(^).^(0)) = ^T(x;L)e'^^ (3.11) 
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where T(x;L) is the sum over all paths of length L that connect to x weighted by 
eicp{—g^Wr} and (■) is the inner product defined on the A^-component internal space. 
This correlation function is not gauge invariant: to make it gauge invariant we could 
add a Wilson line to the correlation function. The arbitrariness of the path we choose to 
connect to x will lead to an arbitrary value of the writhe since the expression in ( p.3| ) 
is manifestly not additive — it depends on the whole curve. 

For a closed polymer loop, however, we can consider the correlation function 

0/3 



^(^"'t(x)-^(0)V^t(o).^(x)) (3.12) 



where a and (3 are the indices of the component field and (■) is the inner product 
defined on the iV-component internal space. This correlation function is gauge invariant 
and, in the N ^ Q limit may be written as 

G^(x)=5:T^(x;L)e'^^ (3.13) 

L 

where now T*" (x; L) is the sum over all closed loops of length L which pass through 
and X weighted by exp{— gf^Wt}. Upon normalizing T*^(x; L) for fixed L we may find 
the probability of a closed, self-avoiding walk passing through and x with a chemical 
potential g"^ for its writhe. For a monomer chemical potential /i close to its critical value 
/ic, G""(x, /i) will scale as [|x| (/i — /ic)'^], where v is the usual correlation length exponent 
of critical phenomena. Inverting the sum in (|3.11|) to find T£(x) will lead to a scaling 
form for the probability 

P^(x;L)(xS(^M^ (3.14) 

and so the radius of gyration of the closed self-avoiding walk will still scale as . In the 
next section we will argue that at the self-avoiding walk fixed point {a.k.a the Wilson- 
Fisher fixed point) the Chern-Simons gauge field is irrelevant, and thus the presence of 
writhe does not change the radius of gyration of a closed polymer loop for g<l. 



4 Renormalization Group Results 

In this section we perform a renormalization group analysis of the complex 0{N) spin 
model coupled to a Chern-Simons gauge field in the — ^ limit. We will show that the 
presence of the Chern-Simons term does not change the long-distance behavior of the 
four point correlation function in (|3.12|) . 
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Our analysis is based on a perturbation expansion around the Wilson-Fisher fixed 
point of the ungauged model. At this fixed point we have the following symbolic recursion 
relations which define the critical behavior for self-avoiding random walks: 

de 



dl 
dn 

H 
du 

H 



e (-?7 + /e(e,M)) 
u fu{e,u). 



(4.1) 



The functions /x can be computed, for instance, by a d = 4 — e expansion The 
scaling relation given in (|4.1| ) for the coefficient e defines an anomalous dimension for 
the scalar fields: if this coefficient is chosen not to rescale, then its exponent is absorbed 
into the ip and ip^ fields. This implies that the scalar propagator no longer scales with its 
naive dimension: at the Wilson-Fisher fixed point the scalar propagator is proportional 
to A;~^+''. This makes the infrared power counting at the Wilson-Fisher fixed point softer 
than that at the Gaussian fixed point where (ip'^ip) ~ kr"^. As we shall see, this difference 



will prove important by keeping the corrections to ([4.1|) from diverging in the infrared. 

In what follows we will consider how the addition of the Chern-Simons gauge 
field modifies these recursion relations for the self-avoiding walk field theory in three 
dimensions. 



4.1 Non-Renormalization of the Gauge Field 

We first show that the gauge propagator is not renormalized. This may not be surprising 
since the gauge field can only couple to closed scalar loops. These scalar loops generate 
diagrams proportional to which vanish in the A^ limit (fig. |l]). To prove that the 
gauge field does not get renormalized, we introduce an auxiliary Hubbard-Stratonovich 
field X through a new term that we will add to the action in ( |3.1CI|) : 



5v 



dh + ^v^l'?/'! 



(4.2) 



Since the x fi^ld enters quadratically into the action, the new partition function will 
simply be proportional to the old one and the equations of motion for ip, and A will 
be left unchanged: 



[(iA][(ix] dip dip^ 



exp 



d^x 



{d — iA) ip — [ji — ixV^uj if 



X 



(4.3) 
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Figure 1: A correction to the gauge propagator with a closed scalar loop which vanishes 
in the iV — > limit. 

Since the scalar field now appears quadratically, we can integrate it out, leaving a 
determinant: 



(4.4) 



J [dA] [dx] det-^ [e {d - iAf - (/i - ix^) 

xexp{-/d%(jLA.VxA+^(V.A)' + ^)) 



In the — > limit, the determinant does not contribute to the A — x partition function. 
The Chern-Simons field A is therefore a free field which keeps its naive dimension in the 
renormalization process. Moreover gauge invariance implies that A must scale as d and 
so, in three dimensions, g"^ does not rescale: 

Though the gauge coupling does not get renormalized, its presence may still renormalize 
the couplings e, /i, and u. We will explore this possibility in the next subsection. 

Superficially it might appear that the same argument presented above for the A field 
could be applied to the scalar field ip. This would lead to the incorrect conclusion that 
ip also keeps its naive dimension upon renormalization. Of course the ip field does get 
renormalized - it acquires an anomalous dimension rj as indicated in ( [4.1|) . The argument 
put forward above cannot be applied to the scalar field as this field's couplings are already 
0{N). The determinant that we were previously able to disregard generates just such 
terms, as can be seen from an expansion of det^[M] = exp(A^tr InM), where M is the 
operator in 



4.2 Perturbation Around the Wilson-Fisher Fixed Point 

We start our renormalization group analysis with the recursion relations for e, n, and 
u at the Wilson-Fisher fixed point of the ungauged 0{N) complex scalar theory in the 
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(b) 






Figure 2: Corrections to the cubic and quartic gauge vertices. 



N ^ limit ( |4.1| ). We will perform a one-loop perturbative expansion in the gauge field 
propagator. 

At zeroth order in our perturbation analysis there is no gauge propagator. The lack 
of a dynamical gauge field severely limits the number of graphs that can modify the 
renormalization group equations: to one loop, there are only three of them (fig. |^). The 
first of these (fig. ||(a)) could contribute a correction to the single gauge field vertex: 

r5it, = / ^^^5^(ki + k, + k3)A.(k3)^^t(k,).^>,)(ki-k,).. (4.6) 

The second and third graphs (figs, ^(b) and 0(c)) could modify the quartic coupling, 

(0) _ f d'^k^ , 

^^^^'^ - J i2nY i2nY (27r)'^ i2nY ^ + ^2 + k, + k,) 

xA,(k3)A(k4)V^^(k2)-^(ki). (4.7) 

In a gauge invariant regularization scheme such as dimensional regularization, fig. |^(a) is 
identically zero. By gauge invariance, the other two graphs must then cancel. Therefore 
to one loop order, the presence of a non-dynamical gauge field does not modify the 
recursion relations given in (|4.1| ) and so leaves the critical behavior of our model 
unchanged. This confirms that A scales as d when A is not dynamical. 



Figure 3: Wavefunction renormalization of the ip field. 
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We now consider terms that are dependent on the gauge propagator. At one loop, 

there is just one new graph in the Landau gauge that could contribute to the scalar 
propagator (fig. ^): 

r^V, = 2//|^0^-(k. + kJv;t(kJ.^-(kJ (4.8) 

x/^(2k.-q).^(-2k.-q),, ' 



(ztt)" q'^ |ki — ^ 

Naive power counting with rj = would have this graph diverging logarithmically in 
three dimensions; however the anomalous dimension of the scalar field [r] ^ 0.016) makes 
this graph finite in the infrared. In fact, r^\^^ vanishes identically because the integrand 
in ( [4. 81) is a contraction of symmetric and antisymmetric tensors. We conclude that the 
presence of the gauge field does not modify the recursion relations for either e or /x to 
one loop. 

We have already seen that there are no corrections to the gauge field vertices at 
one loop order in the absence of the Chern-Simons potential. We now verify that these 
vertices remain unchanged in the presence of the dynamical propagator. Consider the 
one loop corrections to T^^^^. The one loop graphs which renormalize this vertex are 
shown in fig. ^. As a concrete example, we consider the first of these (fig. §(a)) which is 

= 2/ / ^^^S'C-^ + + k3)4(k3)^->(k,).*lk.) (4.9) 



{27tY {27tY 

/■ c^'^q .1 , o IN (2ki + q)i ieijkQj (-2k2 + q)fc 
J [ZttY |ki — qp ^ {k^ + q]"' 



The relevant part of this graph is the part proportional to ki — k2. Terms with higher 
(a) (b) (c) 





Figure 4: Corrections to the single gauge field vertex (H]^) 
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(b) 




(d) 






Figure 5: Corrections to the quartic gauge vertex ( |4.7| ). 

powers of the external momenta generate irrelevant operators and terms with higher 
powers of q generate infrared-finite corrections which do not alter the scaling behavior. 
Leaving these uninteresting pieces behind, we now have 

ri*!* = -io' j '^.^.^/'^^^ + ^^ + ^'')M^^)iH^.)4{^^^) (4.10) 



X 7^(ki-k: 



Once again, a non-zero 7] saves us from having a potentially logarithmic infrared 
divergence for ci = 3. But again, the relevant part of r^^|^ vanishes due to the 
antisymmetry of the Chern-Simons propagator. By power counting, the other two graphs 
(figs. |^(b) and|^(b)) can be shown to generate (at worst) finite corrections. These graphs 
in fact vanish as they are both integrals of odd functions. 

The possible corrections to the quartic gauge coupling ( ^Tzf ) are depicted in fig. |^. 
None of these graphs will contribute to the renormalization of the vertex however: by 
power counting arguments they are all finite at the Wilson-Fisher fixed point. Indeed the 
relevant parts of all of these graphs vanish in Landau gauge due to the antisymmetric 
gauge propagator. 

The one loop results just obtained are natural consequences of gauge invariance. 



6-2?? 
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(a) 




Figure 6: Vanishing corrections to the u vertex. 

We must now consider the u vertex which is marginally relevant at the Wilson- 
Fisher fixed point. We will verify that the fixed point value for u does not lose its 
stability with the introduction of the gauge propagator. To see this, we consider 
the one loop graphs depicted in figs. ^ and |^. Once again power counting tells 
us that none of these graphs can generate a relevant correction to scaling (i.e. an 
infrared- divergent graph). Moreover, the relevant parts of all the graphs in fig. |^ 
vanish, once again due to the antisymmetry of the gauge propagator. The graph 
in fig. 1^ on the other hand, does not vanish and generates a finite correction to u: 

q'i-^dq {d-l)f{d) 



(27r)^ q 



2 



d=3 



«o-^. (4.11) 



P 



Here f{d) is the area of a unit ci-sphere and the persistence length of the chain Lp 
defines the short-distance, ultraviolet cutoff A = L~^. This term decreases the amount 
of self-avoidance in our random walk by a finite amount. Since a finite correction does 
not change the recursion relations in ( |4.1| ), this correction will not change the scaling 
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Figure 7: A finite one loop correction to the u vertex. 



exponent for tlie radius of gyration. Nonetheless, its presence can be felt if the chemical 
potential for writhe g"^ exceeds a certain threshold. When Ag'^L'^Tc^'^ > uq (where Uq is 
the bare repulsion) the repulsive term becomes attractive and the chain collapses as it 
would in a bad solvent. In the case of excluded volume interactions Uq ~ Lp (i.e. fc^T 
per persistence length) and the chain will collapse for (7^>('7r/2)^. 

We have thus shown that the presence of a Chern- Simons field does not change the 
critical behavior of the 0{N) complex scalar model in the — > limit at one loop. 
Using the mappings of sections 2 and 3, we can interpret this results in terms of closed 
self-avoiding walks with a chemical potential for link. The presence of such a chemical 
potential does not affect the statistics of the random walk unless it exceeds a certain 
threshold, at which point the chain collapses. 



5 Average Writhe of a Closed Loop 

We have found that the presence of a small chemical potential for writhe does not change 
the scaling behavior of a self-avoiding random walk. Here we study the first two moments 
of the writhe, (Wr) and (Wr^). 

In the absence of a chemical potential g"^, a rigorous result establishes that 
the ensemble average of the magnitude of the writhe must scale at least as fast 
as the square root of the chain length [0. Numerical results saturate this 
bound: (|Wr|) ~ |T^. We therefore expect the average of the writhe 



squared to scale with the length of the polymer. A simple estimate of {Wr^) 
can be obtained by assuming that the tangents are perfectly correlated within a 



persistence length and completely uncorrelated at larger distances We have 
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1 A.oi„-F / f dR{s,)^dR{s,)\ R{s^)-R{s,] 



(Wi^) = „^ [[dR]e ^ [ dsids2ds3ds4 

lOTT^Z J J 



dsi ds2 I |R.(si) — R(S2) 



3 



f dRjs,) dRjsM R(s3)-R(g4) , . 

I dSs dS4 / |R(S3) -R(S4)|3 ^ ^ 



where the free energy cost is just a bending stiffness: 



F[R] = - ds -—]. (5.2) 



2 J \ds ^ 

Our simphfying assumption is 

/ dRi{s) dRj{s') \ _ j Sij for |s - s'l < Lp 

\ ds ds' / ] otherwise 



(5.3) 



Using this tangent-tangent correlation function and neglecting the correlations between 
R and dR/ds, we find 



[Wr^) ~ J dsids2dszdsi 



L'i / dsTds2-, — 7- , 
Py '|R(si)-R(s2 

Lp / dSids2-pr. r^- 

Li\Si-S2\^ 



(R(si)-R(g2)) (R(s3)-R(g4)) 
|R(S1)-R(S2)|3'|R(S3)-R(S4)|3 

1 



~ Hi; - 1) ~ i;- <^-^) 

Here we have used the statistics for a Gaussian chain in the last step: (|R(s) — R(s')p) ~ 
Lp|s — s'\. Using the self-avoiding scaling behavior (|R(s) — R(s')p) ~ L'^'^'^\s — s'\'^'^ 
gives precisely the same result. We thus see that short-distance correlations saturate the 
rigorous lower bound. 

Our long wavelength theory makes a different prediction for the scaling behavior of 
(Wr^). In section ^, we found the correlation function G"^(x;/i) for a closed polymer 
loop with a chemical potential /i per monomer ( [3.131 ). Using this function as a starting 
point, we determine the scaling behavior of the moments of the writhe. After taking the 
continuum limit for large L, we recover the fixed length generating function using an 
inverse Laplace transform: 

T^(L) = -1. /G''(/i;/)e-^(^-^)d/. (5.5) 
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Figure 8: The writhe is dominated by short plectonemic regions separated by distances 
larger than the persistence length. 



where the contour of integration encloses all the poles of the integrand and 

1 



G {iJ,;g 



lim-(^t(o).^(0; 



(5.6) 



is the generating function for all closed loops. From this function, the second moment of 
the writhe for the closed paths is easily obtained: 



dig 



lnT^(x;L) 



2^2 



(5.7) 



92=0 



Explicit evaluation of the integral in (|5.5|) is difficult, but fortunately we are only 
interested in the way T*-^ scales with L. We assume a scaling relation for the correlation 
function, 

(5.8) 



''^-L--'G''{f,L;g') 



where a is the specific heat exponent and we have used the fact that g"^ does not rescale. 



We can now obtain the scaling behavior of T : 



(5.9) 



As we saw in the last section, the chemical potential for writhe g"^ does not get 
renormalized. What is more, this chemical potential does not appear at this order in 



the scaling exponent of T . It follows that all of the L dependence drops out of ( [5.7| ) 
when the derivatives of the logarithm are taken. We conclude that (Wr^) does not scale 
linearly with the length of the chain; at best, the length-dependence enters as some 
power of a logarithm. Another paradoxical result is obtained when we compute the 



15 



average writhe in the presence of a chemical potentiah 

(Wr) = -:j7^1n T^(L) ~ (5.10) 

The length dependence drops out just as it did for the second moment. A chemical 
potential for writhe does not seem to require the total amount of writhe to grow with 
the chain length. 

These results, though at first surprising, are consistent with existing results. By 
construction, our field theory only applies to the behavior of the random walk on scales 
longer than the persistence length. We have seen that the writhe is concentrated on 
shorter length scales, and so we would not expect to see it in this field theoretic treatment. 
A picture emerges of short plectonemic regions separated from each other by at least the 
polymer persistence length (see fig. |). These plectonemic regions each contribute to the 
writhe of the polymer, but these contributions appear only as matching conditions in our 
long wavelength theory. We would therefore expect that in the presence of a chemical 
potential for writhe, the average writhe of the chain will scale with the chain length: 

(Wr) ^ ^(Wr)pi, (5.11) 

where {Wr)pi is the average writhe per plectonemic segment. Likewise, the second 
moment is dominated by the small plectonemic regions: 

{Wr") ^ ^{Wr\,. (5.12) 

It might seem strange to talk about the writhe of a short segment of the polymer when 
writhe is a global property of the entire chain. But as can be seen from the | R(s) — R(s') | ~^ 
factor in ( |2.3|) , the writhe is dominated by the correlations between regions that are close 
together in space. 

6 Conclusion 

In this paper we have investigated the topology of self-avoiding random walks. Our 
goal was to understand how the linking number constraint changes the conformational 
statistics of closed polymers. After integrating out the twist degree of freedom, this 
problem was shown to be equivalent to that of studying the statistics of a polymer with 
a chemical potential for writhe. By mapping this modified problem onto the ^ limit 
of an 0{N) complex scalar field theory with a Chern-Simons gauge field, we were able 
to find that the presence of a small chemical potential for writhe does not change the 
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behavior of our polymer in the scahng hmit: that is, the renormahzation group equations 
for the ungauged scalar field model that describes the usual self-avoiding walk do not get 
modified by the introduction of the Chern-Simons term. We also saw that a sufficiently 
large chemical potential, g'^>{7i/2), collapses the polymer. 

We carried out our perturbation expansion to one loop. The anomalous wavefunction 
renormahzation of ip at the Wilson-Fisher fixed point kept all of the possible corrections 
to scaling finite; this result should be compared to that of naive power counting where 
logarithmic divergences occur. In view of this, it is reasonable to believe that the 
renormahzation group equations do not change to all orders in the loop expansion. It 
would be interesting to carry out a careful analysis of the power counting at the Wilson- 
Fisher fixed point. It may also be possible to adopt the wealth of knowledge on anyons 
2^ to study the critical behavior of our model as well as the use of supersymmetry 



to study self-avoiding walks. 

Additionally, we investigated the scaling behavior of the writhe itself. We found that 
the long wavelength contribution to the writhe does not scale with its length. The writhe 
escapes the notice of the scaling theory by hiding at small scales. We find our results to 
be consistent with a picture whereby the writhe is confined mainly to short plectonemic 
segments separated by at least the bend persistence length. 

Our overall results indicate that the effects of topology do not change the general 
scaling properties of polymers: constrained random walks belong to the same universality 
class as unconstrained ones. The effects of topology seem to enter only on small scales 
through the chemical and physical details of the chain. 



Recently the winding number statistics of directed polymers have been analyzed 
by adding an auxiliary 2 + 1 dimensional gauge field. We note that the Chern-Simons 
formalism presented here can be used to study the linking number statistics of an 
ensemble of closed, isotropic, polymer loops, as can be found, for instance, in the 
kinetoplast DNA of trypanasomes |2^. The rheological properties of these so-called 
"olympic gels" are intimately connected to their topology. Thus, via bulk experiments, 
the linking number statistics could be explored. 
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